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Abstract 

Monte Carlo Renormalization Group (MCRG) methods were designed to study the non- 
perturbative phase structure and critical behavior of statistical systems and quantum field theories. 
I adopt the 2-lattice matching method used extensively in the 1980's and show how it can be used 
to predict the existence of non-perturbative fixed points and their related critical exponents in 
many flavor SU(3) gauge theories. This work serves to test the method and I study relatively well 
understood systems: the Nj = 0, 4 and 16 flavor models. The pure gauge and Nf = 4 systems are 
confining and chirally broken and the MCRG method can predict their bare step scaling functions. 
Results for the A'^^ = 16 model indicate the existence of an infrared fixed point with nearly marginal 
gauge coupling. I present preliminary results for the scaling dimension of the mass at this new 
fixed point. 
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I. INTRODUCTION 



The Monte Carlo Renormalization Group (MCRG) methods, based on Wilson's renor- 

980 's to study the 
The 2-lattice 
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malization group theory, were developed and used extensively in the_ 
critical properties of spin and gauge models 
matching MCRG proved to be particularly useful to calculate the (3 function of asymptot- 
ically free theories, like quenched QCD jj, 0, Ig]. The approach has all but been forgotten 
in the last 20 years as lattice QCD calculations focused on spectral and other experimen- 



tally measurable quantities 
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coupled beyond-Standard Model physics 
the issues and recent lattice results. A basic discussion of the physica 
remarks on expectations for lattice simulations were presented in Ref. 



picture with some 



follow the Wilson renormalization group RG language used in 



3l|. In this paper I 



3l|. 



SU(3) gauge models with Nf fermions in the fundamental representation can have very 



different phase structure depending on the number of fermions [26|. If Nf > 16 asymptotic 
freedom is lost, the gauge coupling is irrelevant and the continuum theory is free. For small 
Nf the g = 0, m = Gaussian fixed point (GFP) is the only critical fixed point (FP). The 
theory is asymptotically free, confining and chirally broken. Somewhere around Nf ^ 10 the 



gauge coupling develops a new FP at g* [29|,l32]. At g* the gauge coupling is irrelevant, it 



is an infrared FP (IRFP). The continuum limit defined in the basin of attraction of this IRFP 
is neither confining not chirally broken; it is conformal when m = 0. This conformal phase is 
expected to exits all the way to Nf = 16. Identifying the lower end of the conformal window 
and the critical properties of the IRFP are the main issues of recent lattice simulations. The 
MCRG method was designed to answer these kind of questions and in this paper I present 
the first such study in Nf = 4 and Nf = 16 flavor SU(3) theories. I also investigate the 
pure gauge SU(3) model where it is possible to do high statistics, large volume simulations. 
I have chosen these models as the expected phase structure is rather well known, so I can 
use them to calibrate and test the method. My eventual goal is to extend these studies to 
other flavor numbers or fermions in different representations. 

Since MCRG has been used very little in the last 20 years, I devote Sect JIIII to the basic 
description of the 2-lattice matching MCRG. The method allows the determination of a 
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sequence of couplings /3o, ... with lattice spacings that differ by a factor of s between 

consecutive points, a(/5„) = a(/5„_i)/s. s is the scale change of the RG transformation, 
s = 2 in this study. This sequence is analogous to the step scaling function defined in the 



Schrodinger functional (SF) method 33|, |3J, |35|] , but in MCRG it is defined through the bare 



couplings. To emphasize this difference I will use the notation Sfe(/5„; s) = /?„ — for the 
bare step scaling function instead of the more traditional a{u; s) used in the SF approach. 

The sequence Po, Pi, ... can be used to determine the renormalized running coupling 
in theories that are governed by the GFP if at the weak coupling end of the chain a renor- 
malized coupling, like the SF g^, is calculated and connected to a continuum regularization 
scheme, while at the strong coupling end some physical quantity is used to determine the lat- 
tice scale. I do not pursue this calculation here, though I will compare results for Sh{P', s = 2) 
from SF and MCRG in Sect. HVAl 

The rest of the paper is organized as follows. Sect. [TTl summarizes the perturbative 
picture of these many fermion theories. Sect. IIIII describes the 2-lattice matching method 
and defines the RG block transformation used in this work. The numerical simulations and 
results are discussed in Sect. IIVI The technical aspects of MCRG are described in detail 
for the pure gauge SU(3) theory in Sect II V A] as it serves to justify the approach used for 
the Nf = 4 and Nf = 16 models in Sects. IIVBI and llVCi I use nHYP smeared staggered 
fermions in this work. I present some basic properties of 4 flavor staggered fermions, together 
with the MCRG calculation of the bare step scaling function in Sect. IIVBI The existence of 
an IRFP requires the re-evaluation of the MCRG method. This, together with preliminary 
results for the scaling exponent of the mass in the Nf = 16 model are presented in Sect. 
HVCl 

II. THE PERTURBATIVE PICTURE 

Before discussing the MCRG method and the numerical results I briefly summarize the 
perturbative picture. The universal 2-loop /3 function for SU(3) gauge with Nf fermions in 
the fundamental representation is 



d\og{jj'^) 167r2 (167r 



2^2• 
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b2 = -102 + ^iV/. 

For Nf < 16.5 the 1-loop coefficient bi is negative, the gauge couphng is relevant at 
the g = Gaussian FP, the theory is asymptotically free. Dimensional transmutation is 
responsible for mass generation. The energy scale changes by a factor of 2 between couplings 
gi and g2 if 

At one loop level this leads to a constant shift in /? = 6/(7^ and the bare step scaling function 
is 

s,{P,;s = 2) = P,-P2 = -^^b, (1-loop). (3) 

For small fermion numbers the higher order terms are small, the j3 function is expected 
to remain negative. Lattice simulations indicate that for Nf < 8 the system is confining and 
chiral symmetry is spontaneously broken. For Nf > 8 the 2-loop (3 function develops a zero 



at g* ^ Banks-Zaks FP 26!]. At this new FP g is irrelevant, it is an IRFP for the gauge 
coupling. The infinite cut-off limit in the vicinity of g* is conformal. 

When the perturbatively predicted g* is large, higher order or non-perturbative effects 
can destroy the existence of the IRFP. Analytical considerations and numerical simulations 



suggest that the bottom of the conformal window is around Nf ^ 10 [29|, |32l]. At Nf = 16, 
the largest flavor number that is still asymptotically free, the Banks-Zaks FP occours at 
a small value g* ~ e, perturbation theory could correctly describe the conformal phase. 
At the IRFP there is only one relevant operator, the mass. Its scaling dimension (critical 
exponent) is close to its engineering one ~ 1 + 0{e), while the scaling dimension of the 
gauge coupling is i/g ~ — e^. The slope of the /3 function at g* predicts the exponent 

Pig') = -y^g' - g*') + o^g' - g*^) , (4) 

^ ' (5) 



'2 



Eq. [2] now gives 



gf - g*' = [gl - g*')2~''^ , (6) 

s,(/?i; s = 2) = A - /52 = (/?2 - /3*)(2-2^« - 1) , (7) 
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if Pi — (3* , — ^ 1- For 16 flavors perturbatively i/g ^ —0.01, the gauge coupling is 
almost marginal. For smaller Nf or higher representation fermions \yg\ can be larger, though 
both numerical and analytical considerations find that \yg\ remains small even at the bottom 



of the conformal window 
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31 



36|. 



The mass is a relevant operator both at the GFP and at the IRFP, with critical value 
m* = 0. Under a scale change s = 2 it changes as 

mi = m22~^™ , (8) 

where l/i/m = is the critical index of the mass. 

III. THE MCRG METHOD 

The Wilson RG description of statistical systems is a very effective approach to describe 
the phase diagram, calculate critical indices, and in case of lattice discretized quantum field 
theories, understand the infinite cut-off continuum limit of these models. There are many 
books and review articles written about the subject. I do not attempt to explain Wilson RG 



here, I only summarize _t 
of this paper are Refs 



uhe main points. Two reviews that could be useful for other parts 

"0. 

In the inherently non-perturbative Wilson RG approach one considers the evolution of 
all the possible couplings under an RG transformation that preserves the internal symme- 
tries of the system but integrates out the cut-off level UV modes. The fixed points of the 
transformation are characterized by the number of relevant operators, i.e. couplings with 
positive scaling dimensions that flow away from the FP. Irrelevant couplings have negative 
scaling dimensions and they flow towards the FP. The IR values of irrelevant operators are 
independent of their UV values. Continuum (or infinite cut-off) limits can be defined by 
tuning the relevant couplings towards the FP, thus controlling their IR value. The number 
of relevant operators and their speed along the RG flow lines are universal, related to the 
infrared properties of the underlying continuum limit. On the other hand the location of 
the FP is not physical, in fact different RG transformations have different fixed points. 

In quantum field theories the best understood fixed points are at vanishing couplings 
(Gaussian FPs) as they can be treated perturbatively. For example the GFP of the 4 
dimensional SU(3) pure gauge model has one relevant operator, the gauge coupling, and no 
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other FP of the model is known to exist. The Gaussian FP of 2-flavor QCD has two relevant 
operators, the mass and the gauge coupling. Gauge theories with many flavors can develop, 
in addition to the GFP, a new flxed point where only the mass is relevant (Banks-Zaks 
infrared flxed point) 26|. These new FPs are rarely in the perturbative region and to study 
their existence and properties is the main motivation for this paper. 



A. The 2-lattice matching MCRG method 

Consider a d-dimensional lattice model with action S{Ki). {Ki} denotes the set of all 
possible couplings, though in a typical lattice simulation only a few of them are non-zero. 
The system is characterized by one or more length scales, like the correlation length ^, 
inverse quark masses, etc. In numerical simulations we always deal with flnite volume and 
for now I assume a hypercubic geometry with linear size L = La. The flrst step of a real 
space renormalization group block transformation is to deflne block variables. These new 
variables are deflned as some kind of local average of the original lattice variables and for 
a scale s > 1 transformation they live on an L/s lattice. By integrating out the original 
variables while keeping the block variables flxed one removes the ultraviolet fluctuations 
below the length scale sa. The action that describes the dynamics of the block variables 
is usually much more complicated than the original one, but if s is much smaller than the 
lattice correlation length, the long distance infrared properties of the system are unchanged. 
After repeated block transformation steps the blocked actions describe a flow line in the 
multi-dimensional action space 

{K,} ^ {irf {irf } ^ .... , (9) 

where {K^^^} denotes the couplings after n blocking steps. While the physical correlation 
length is unchanged, the lattice correlation length after n blocking steps is 

^(n) _ s-'^^(O). (10) 

The RG can have flxed points only when ^ = oo (critical) or ^ = (trivial). We are, 
of course, interested in the former one. Near the critical flxed point the linearized RG 
transformation predicts the scaling operators and their corresponding scaling dimensions. 

It is easy to visualize the renormalization group flow lines when there is only one relevant 
coupling at the flxed point, as illustrated in Fig. [H The sketch depicts the flow lines in 
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K K' Kq 

FIG. 1: Sketch of the RG flow around a FP with one relevant operator. The coupling pair (K, K') 
indicates matched couplings whose correlation length differ by a factor of s. 

the parameter space {Kq^Ki}, where for simplicity I assume that the critical surface is at 
Kq = 0. Flow lines starting near the critical surface approach the fixed point in the irrelevant 
direction(s) but flow away in the relevant one. After a few RG steps the irrelevant operators 
die out and the flow follows the unique renormalized trajectory (RT), independent of the 
original couplings. If we can identify two sets of couplings, {Ki} and {K-}, that end up 
at the same point along the RT after repeated blocking steps, we can conclude that their 
correlation lengths are identical. If they end up at the same point along the RT but one 
requires one less blocking steps to do so, according to Eq. [10] their lattice correlation lengths 
differ by a factor of s. This is also illustrated in Fig. [1] From K' one needs 3 while from 
K one needs 4 RG steps to reach the same point of the RT (up to small corrections in the 
irrelevant direction(s)), therefore ^' = ^/s. This gives the bare step scaling function with 
scale change s. The two lattice matching 0, is a numerical method to identify {K, K') 
pairs. 

In order to identify a pair of couplings {K, K') with ^' = ^/s we have to show that after 
n and {n — 1) blocking steps their actions are identical, S{K^"'') = S{K'^^ ^''). It is quite 
difficult to calculate the blocked action, but fortunately we do not need to know the actions 
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explicitly to shows that they are identical. It is sufficient to show that the expectation 
values of every operator measured on configurations generated with one or the other action 
are identical. Furthermore it is possible to create a configuration ensemble with Boltzman 
weight of an RG blocked action by generating an ensemble with the original action and 
blocking the configurations themselves [1]. This suggests the following procedure for the 
2-lattice matching: 

1. Generate a configuration ensemble of size L'^ with action S{K). Block each configu- 
ration n times and measure a set of expectation values on the resulting [L/ s'^Y set. 

2. Generate configurations of size {L/sY with action S{K'), where K' is a trial coupling. 
Block each configuration n — 1 times and measure the same expectation values on the 
resulting {L/s"'Y set. Compare the results with that obtained in step 1. and tune the 
coupling K' such that the expectation values agree. 

A few basic comments are in order: 

a) Since we always compare measurements on the same lattice size, the finite volume 
corrections are minimal and even very small lattices can be used. 

b) It is not necessary to work on lattices that are larger than the correlation length of the 
system, nor does it matter if we are in the confined or deconfined phase of the system. 

c) If the fiow lines follow the unique RT, even one operator expectation value is sufficient 
to find the matching coupling, all other operators should give the same prediction. In 
practice we can do only a few blocking steps and the fiow lines might not reach the 
RT. That will be refiected by different operators predicting different matching values. 
The spread of these predictions measure the goodness of the matching. Increasing 
the number of blocking steps improves the matching, and when the RT is reached, 
consequetive blocking steps predict the same matching couplings. 

d) The location of the fixed point and its renormalized trajectory in the irrelevant di- 
rections depend on the block transformation. Block transformations that have free 
parameters can be optimized so their RT is reached fast and the matching is reliable 
after a ver y fe w RG steps. This optimization proved essential in previous applications 
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e) Since we can match by comparing local operators, the statistical accuracy is usually 
acceptable even with small configuration sets. 

If the FP has two relevant operators, the matching proceeds similarly but one has to tune 
2 operators. In practice this is much more difficult than the tuning of a simple coupling. It 
is frequently easier to fix one of the relevant couplings to its FP value and proceed with the 
matching in the second relevant coupling as described above. 

I will illustrate the above points in Sect II VI 



B. The renormalization group block transformation 

I chose a scale s = 2 block transformation, similar to what was used in Refs. 2|, |5j 
where Proj indicates projection to SU{3). The parameter a is arbit^ar^ and can be used 



roitrary 



to optimize the blocking. The block transformation used in Refs. p, |6|, |37| had a fixed, 
1 — a = a/Q, but instead of projecting to SU(3) the blocked link was allowed to fiuctuate 
around V^n.^i, depending on a free parameter. In my experience the two block transformations 
are very similar. 

In principle one can define an RG transformation for fermions as well. However it is 
easier to do the RG transformation after the fermions are integrated out, i.e. when the 
action depends on the gauge fields only. 

The role of the parameter a is to optimize the block transformation. While the critical 
surface of a system is well defined, the location of the fixed point itself is not physical, it can 
be changed by changing the RG transformation. It is important to optimize the blocking so 
its FP and RT can be reached in a few steps. The optimal blocking is characterized by 

1. Consistent matching between the different operators: along the RT all expectation 
values should agree on the matched configuration sets. Any deviation is a measure 
that the RT has not been reached. 

2. Consecutive blocking steps should give the same matching coupling. When they pre- 
dict different values, one can try to extrapolate to the FP using the first non-leading 
critical exponent. 



In the next Section I will show that both of the above conditions can be satisfied in 
numerical simulations if the blocking parameter is optimized. 



IV. SIMULATIONS 



A. SU(3) pure gauge theory 



At the Gaussian g = FP of the pure gauge SU(3) model the gauge coupling is relevant, 
the theory is asymptotically free. According to Eq. [3] 1-loop perturbation theory predicts 
that the bare step scaling function is constant, independent of the gauge coupling. The 2- 
loop corrections are small in a wide range of coupling, sf^'^'^*'* fa 0.59 is a good approximation. 



The bare step scaling function was studied in Refs. [4, l5|, la, Il0[ with the 2-lattice matching 
method. Here I repeat some of those calculations with a different block transformation and 
extend them to larger volumes and statistics. Where they overlap, the results I present 
below are consistent with the original calculations. This section mainly serves as a test of 
the method. 

I generated 200-300 independent configurations at several coupling values with the Wilson 
plaquette gauge action and calculated the bare step scaling function Sb{l3; s = 2) matching 
32^ volumes on to 16'^ , and also 16^ volumes to 8^. The 32'^ volume can be blocked up to 
4 times and compared to the 16^ volume that is blocked up to 3 times. At each blocking 
level I measured 5 operators: the plaquette, the 3 6-link loops and a randomly chosen 8-link 
loop. 

Figure [2] illustrates the 2-lattice MCRG. The plot shows the matching of the plaquette 
with the s = 2 renormalization group transformation of Eq. [11] and blocking parameter 
a = 0.65. The 32^ volume simulations were done at /5 = 7.0, and the cyan, blue and red 
symbols are the values of the blocked plaquette after 2, 3 and 4 blocking steps. The solid 
curves interpolate the plaquette values, measured at many couplings on 16* volumes, after 
1, 2 and 3 blocking steps. The 32^ data match the 16^ values at jS' = 6.49 for all blocking 
levels. The final blocked volume is 2^, but finite size effects are minimal as one always 
compares observables on the same volume. 

The matching can be repeated with different operators and RG transformations. Figure 
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6 7 8 

FIG. 2: The matching of the plaquette for pure gauge SU(3) theory. The simulations were done 
on 32^ volumes at /? = 7.0 (symbols) and 16^ volumes at many coupling values (solid interpolating 
lines). The configurations were blocked with s = 2, a = 0.65 parameter block transformation 2(1) 
(cyan), 3(2) (blue) and 4(3) times (red). 

El shows the difference between the matched couplings, 

A(3 = (3-i3', (12) 

as the function of the blocking parameter a for the 5 different operators at the last 3 blocking 
levels for (3 = 7.0. The plots show two trends. First, the spread of the predicted A/5 values 
from the different operators decrease with increasing blocking level, signaling that the RG 
flow lines are approaching the RT of the block transformation. Second, the dependence of 
Af3 on the blocking parameter decreases with increasing blocking levels suggesting a unique 
value for A/5 in the rii, —>■ oo limit. 

Figure S] summarizes the plots of Figure [31 The average matching A/5 values are plotted 
as the function of a for the last three blocking levels. The "error bars" show the standard 
deviation (spread) of the predicted values, therefore they represent the systematic errors 
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FIG. 3: Matching at /? = 7.0 from 32^ to 16^ lattices. The matching values A/3 as the function of 
the blocking parameter a are shown for the 5 different operators measured. Left panel: blocking 
level rife = 2(1); middle panel n;, = 3(2); right panel n^, = 4(3). 



of the matching procedure. The statistical errors are small, comparable to the systematic 
errors only at the last blocking level at the best matching around a = 0.65. The 3 different 
blocking levels converge around a = 0.65, the same value where the spread from the different 
operators is minimal, predicting the relation between the lattice spacings a{f3' = 6.485) = 
2a(/5 = 7.0). In the Uf, ^ oo limit the quantity A/3(/5) = (3 — (3' is the bare step scaling 
function Sb{P; s = 2), analogue to the step scaling function of the renormalized coupling 
used in the SF formalism. In the following I will use the intersection of the last two blocking 
levels to identify Sft(/?) as it is usually less sensitive to the statistical errors than the spread 
of the individual operators. 

Table [T] summarizes the results at different couplings and volumes, together with the 
optimal blocking parameters a. The data indicate consistency between the different volumes 
and increasing blocking levels. This observation will be important in the study of the Nf = 4 
and 16 systems where matching is done on 16'^ — * 8^ lattices only. The uncertainty of the 
predictions from the 16^ 8^ are considerably larger than from the larger volumes. This is 
not statistical, rather reflects the fact that the systematical errors of the matching after 3(2) 
blocking steps are larger than after one more blocking level. It is possible that a different 



block transformation would give better matching. I have not been ab 
s = 2 transformation to make it better. Adding HYP smearing 



e to modify the scale 
] before constructing 
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FIG. 4: Matching at /3 = 7.0 from 32^ to 16^ lattices. The average matching values A/3 as the 
function of the blocking parameter for the last 3 blocking levels are shown. Note that the "error 
bars" denote the spread of the predictions for the 5 different operators used and thus represent 
systematical errors. 

the blocked links pulls the RT closer, but at the same time reduces the dependence of the 
expectation values on the couplings, thus increasing the errors. It would be worthwhile to 
try combining the block transformation with a simple APE smearing, or use a variation of 
the scale s = \/3 transformation of Refs. 10, lUj that would allow more blocking steps from 
the same lattice size. 

The quantity Sb{P; s) is the bare step scaling function for scale change s = 2. One can 



predict its value using physical observables like the Sommer parameter ro[39| or the critical 
temperature Tc . The SF calculation or the recently proposed new method to calculate the 
renormalized coupling based on Wilson loop matching (WL) [40] can also predict Sf,(/5). 

In Figure [5] I compare the MCRG result for the bare step scaling function with predictions 
from other methods. Note that I show errors only for the MCRG results. I used the 



interpolating formula from Ref. [4l| to find /?') pairs where rg/a differ by a factor of 
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Q^opt 


Sfe, L = 32 
ub = 3(2) 


Sb, L = 32 
ub = 4(3) 


Sb, L = 16 
Ub = 3(2) 


6.0 


0.71 






0.365(4) 


6.2 


0.72 






0.410(7) 


6.4 


0.71 


0.451(12) 


0.448(10) 


0.468(16) 


6.6 


0.69 


0.488(15) 


0.483(5) 


0.496(13) 


6.8 


0.66 






0.511(19) 


7.0 


0.66 


0.517(27) 


0.515(6) 


0.516(10) 


7.2 


0.63 






0.536(26) 


7.4 


0.61 


0.548(38) 


0.571(6) 


0.575(42) 


7.8 


0.60 


0.558(34) 


0.575(5) 


0.573(42) 



TABLE I: The bare step scaling function for the pure gauge SU(3) system. The second cohimn 
hst the optimal blocking parameter. The third and fourth columns are results from simulations 
on 32^ volumes matched to 16'^ after 3(2) and 4(3) blocking steps. The last column shows results 
from 16^ volumes matched to 8'^ after 3(2) blocking steps. 



2, while for Tg I used the Nt = 8 and 4, and Nt = 12 and 6 transition temperatures 
from Ref. |42j|. In case of the SF and WL calculations I attempted to find matching 
(/?, /3') pairs by identifying bare couplings where the renormalized SF couplings are related 
as g'^{L) = (7'^ (21/). In principle this relation should be taken in the L — 00 limit but the 
numerical data do not show significant finite volume effects. The predictions in Figure [5] use 
data from the 1-loop improved SF [34] . The bare couplings used in the 2-loop improved SF 
paper do not match close enough to use them in this analysis 35 1. 

All the above calculations use the Wilson plaquette gauge action, so in the scaling regime 
they should give the same prediction. It is very satisfying to see the agreement between 
MCRG, To and Tc even at relatively strong couplings. In the range (3 G (6.0, 7.0) the 
predicted values differ considerably from the 2-loop perturbative results. It is difficult to 
measure tq or Tc at much finer lattice spacings and show perturbative scaling for them. On 
the other hand both the SF, WL and MCRG methods approach s^^'^^\ but the latter one 



only atf3> 7.0. 
analyzed in Ref. 



he relatively large difference between the SF and tq data was discussed and 



4l| where it was also noted that the 2-loop improved SF shows significantly 
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FIG. 5: The bare step scaling function si,{P;s = 2) for the pure gauge SU(3) system as predicted 
by different methods. The 1-loop perturbative prediction is s^^'^'^*^(s = 2) = 0.59. 



smaller scaling violations relative to tq. 

Based on the results presented in this section the 2-lattice MCRG matching method could 
be competitive with other methods in determining the running coupling of asymptotically 
free theories when it is combined with an independent definition of the renormalized coupling 
in the weak coupling regime. 



B. Nf = 4 flavor model 

The 4 flavor SU(3) gauge theory is expected to be confining and chirally broken even 
at large gauge couplings. At the Gaussian g = 0, m = FP both the mass and the 
gauge coupling are relevant operators, the model is asymptotically free. Perturbation theory 
predicts s^^'''^\s = 2) = 0.45. 



The results 1 present here were obtained using nHYP smeared staggered fermions 43 1. 



1 chose nHYP smearing as it significantly reduces taste breaking of staggered fermions and 
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therefore even in strong coupling has manageable lattice artifacts. 



1. The nHYP staggered action 



Very little is known about the 4-flavor system with nHYP or HYP smeared fermions. 



The finite temperature phase transition of 
partial-global Monte Carlo update in Ref. 



;he HYP smeared model was studied with the 



44| . The phase transition in the chiral limit is 



expected to be first order, most likely extending to finite mass before turning into a crossover 
at large masses. Simulations with thin link staggered fermions have confirmed this, finding 
a strong discontinuity even at fairly large quark masses. The conclusion of Ref. was 
quite different: we found no signal for discontinuity, the phase transition appeared to be a 
crossover both for Nt = 4 and 6 even at fairly small masses. The updating technique used 
in Ref. 



44l | was not efficient enough to pursue much larger volumes, and we did not continue 
our investigation of the Nf = 4 system. 
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FIG. 6: The condensate and the disconnected chiral condensate on 8^ x 4 volumes at m = 0.04 in 
the Nf = 4 theory. 

The nHYP smeared action is nearly identical to the HYP smeared one, but nHYP is 
different iable and the efficient molecular dynamics update can be used with it j43|. I have 



confirmed the raw data of Ref. 



44l | with the nHYP action, and extended it further toward 



the strong coupling region. Figure [6] shows the condensate {ijjijj) and the disconnected chiral 
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FIG. 7: The condensate and the disconnected chiral condensate on 12'^ x 6 volumes at m = 0.05 
in the Nf = 4 theory. 



susceptibility 



45| 



U 



(13) 



on 8^x4 lattices at m = 0.04. The condensate is almost identical to Figure 5 of 4J], a smooth 
function of the gauge coupling with no obvious sign of discontinuity. The disconnected chiral 
susceptibility has a strong peak at /3 = 4.4, signaling a crossover. In contrast, the data with 
thin link staggered fermions at Nt = 4 show a discontinuity in the condensate A{iljip) ^ 0.2. 
Figure[7]is the same as Figure[6]but on 12^ x 6 lattices at m = 0.05. Again, the condensate is 
smooth, the susceptibility suggests a crossover around (3 = 5.0. Obviously much more work 
is needed to determine the transition temperatures and the order of the phase transition 
accurately. Larger spatial volumes might sharpen the transition, but in any case the endpoint 
of the first order line with nHYP fermions occurs at much smaller masses than with the thin 
link action. 

To set the scale I measured the static potential on 16^ lattices used later in the MCRG 
study. I found ro/a = 5.8(3) at f3 = 5.4, m = 0.01. Ref. quotes ro/a = 3.34(7) at 
/3 = 5.2, m = 0.04 and ro/a = 2.2(1) at /3 = 5.0, m = 0.10. 
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2. MCRG matching 
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FIG. 8: Matching of the plaquette in the A'^^ = 4 model. The blue and red symbols show the 
plaquette value after 2 and 3 blocking steps on the 16^ volumes at /3 = 6.0, m = 0.01. The lines 
interpolate the plaquette on the 8'^ volumes at several couplings and m = 0.015 (dashed) and 
m = 0.025 (dotted) after 1 and 2 blocking levels. The black line is the unblocked plaquette on the 
8^ lattices. The dashed and dotted curves are barely distinguishable. 

In principle MCRG matching could be done similarly to the pure gauge SU(3) model, but 
since the fermionic model has 2 relevant operators, the matching requires tuning in both 
the gauge coupling and the mass . This is a considerably harder numerical task than a 
single parameter matching. One can reduce this complication by setting one of the relevant 
couplings to its critical value, since then only the other coupling has to be matched. The 
critical value of the mass is m = 0. Simulations in small volumes are possible even with 
vanishing mass. In addition the dependence of the local observables used in the matching 
is so weak on the mass that a small mass in the simulations is also acceptable. Setting the 
mass to zero or to a small value allows matching in the gauge coupling only and the bare 
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a 

FIG. 9: Matching at /? = 6.0 from 16^ to 'S^ lattices in the A'^^ = 4 model. The matching values A/3 
as the function of the blocking parameter for the last 2 blocking levels are shown. Note that the 
error bars denote the spread of the predictions for the 5 different operators used and thus represent 
systematical errors. 

step scaling function can be calculated in the same way as for the pure gauge system. 

To calculate the step scaling function I considered 16^ — > 8^ matching at several gauge 
couplings between (3 G (5.4,8.0) (see Table [IT|) . All the 16'^ configurations except [3 = 5.4 
are in the deconfined, chirally symmetric phase, but that does not matter for the MCRG 
matching method. I have generated 100-150 configurations on the larger volumes and ~ 300 
configurations on the smaller ones, separated by 10 molecular dynamics trajectories. I have 
used the same 5 operators in the matching as in the pure gauge SU(3) system. 

On the 16^ lattices I chose m = 0.01. If the critical exponent for the mass were its 
engineering dimension, the matching mass on the smaller volume would be m = 0.02. I 
generated 8^ lattices with m = 0.015 and 0.025 to bracket this value and to check for any 
dependence on the mass. The matching for the plaquette is shown in Figure [HI The blue 
and red symbols represent the plaquette value after 2 and 3 blocking steps on the 16^ lattice 




19 



/3l6 




Sb{s = 2) 


5.4 


0.81 


0.388(10) 


5.6 


0.81 


0.330(13) 


5.8 


0.77 


0.335(20) 


6.0 


0.76 


0.303(25) 


6.4 


0.67 


0.400(29) 


6.8 


0.67 


0.365(42) 


7.2 


0.63 


0.404(39) 


8.0 


0.59 


0.470(53) 



TABLE II: The bare step scaling function Sf,(s = 2) for the 4-flavor simulation. The second column 
lists the optimal blocking parameter a. 
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FIG. 10: The bare step scaling function for 4 flavor SU(3) theory. The dashed line indicates the 
1-loop perturbative value. 
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at /3 = 6.0 while the dashed and dotted hnes interpolate the 1 and 2 times blocked plaquette 
values on the 8^ lattices with m = 0.015 and m = 0.025. For completeness I also show the 
unblocked plaquette on the 8^ volumes (black line) though there is no consistent matching 
at this level. The dotted and dashed curves are barely distinguishable. None of the other 
observables show any dependence on the mass even after 3 blocking steps beyond the fairly 
small statistical errors of the simulations, implying that the system indeed can be considered 
critical in the mass. 

The dependence of the matched values on the blocking parameter a is similar to the pure 
gauge case. The analogue of Figure H] is Figure [9l Since for 16"^ 8"^ matching only two 
blocking levels can be used, one has only two sets of predictions, Uf, = 2(1) and Uf, = 3(2). In 
Sect. IIV A[ Table H] I found that one can safely identify the optimal blocking parameter and 
matching value from the intersection of these blocking levels even on 16^ 8^ matching. 

Table [TTl and Figure [10] summarize the simulation results. The step scaling function 
Sb{s = 2) is consistent with asymptotic freedom and approaches the 2-loop perturbative 
prediction at weak couplings. Just like in the pure gauge system, the data in Table [III could 
be used to determine the running coupling of the 4-flavor model. To do so one needs to 
calculate a renormalized coupling aX (3 = 8.0 and connect it through perturbation theory to 
a continuum scheme. The change of the lattice scale between (3 = 5.4 and (3 = 8.0 can be 
determined form the data, while the lattice spacing can be obtained by measuring a physical 
quantity at some strong coupling. Improved block transformation or larger volumes would 
allow a more precise determination of Sb(/?). Higher statistics especially at the larger (3 
values, would also help. 

C. Nf = 16 flavor model 

The 16 flavor SU(3) model is still asymptotically free, but 2-loop perturbation theory 
predicts the emergence of an IRFP at weak gauge coupling Continuum limits can be 
defined both at the Gaussian FP and at the new IRFP. In the former case both the gauge 
coupling and the mass has to be tuned towards the FP, in the latter one the gauge coupling is 
irrelevant, only the mass has to be tuned to m = 0. There is no confinement or spontaneous 
chiral symmetry breaking in the weak coupling phase, the continuum massless theory at the 
IRFP is conformal. The conformal phase was first identified in Ref. 
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It is generally believed that in the strong coupling the lattice model is confining and 
chirally broken, so there has to be a (bulk) transition separating the strong and weak coupling 
phases. Since the bulk transition is a lattice artifact, it is probably not associated with 
critical behavior or continuum quantum field theory. Most likely it is a first order phase 
transition at m = and it might extend to m > before turning into a crossover [si], H] ■ 
One should mention that some numerical results indicate that this confining phase might 



not even exists 



48|. 



1. MCRG around an infrared fixed point 

Assuming that the simulations are done in the conformal phase, the behavior we expect 
from MCRG depends on whether we study the critical (m = 0) or the m 7^ phases: 

• On the m = critical surface at very weak coupling the system is in the attractive 
region of the Gaussian FP. 2-lattice matching could reveal the running of the gauge 
coupling, sf^^^"* = 0.016 from the 1-loop /? function, though numerical simulations 
probably will not be close enough to the Gaussian FP to see this behavior. It is much 
more likely that the RG fiow will be determined by the Banks-Zaks IRFP. At this FP 
all operators are irrelevant when m = 0. The 2-lattice matching is designed to map 
out the fiow speed in the relevant direction, so we have to re-evaluate the method 
when there is no relevant operator. 

When all operators are irrelevant, they all fiow into the FP according to their scaling 
dimensions. In the rih 00 limit all expectation values approach their FP value 
independent of the bare couplings, so matching is meaningless. At finite rib the RG fiow 
can pick up the " least irrelevant" operator. If there is one operator with a nearly zero 
scaling dimension matching can make sense: after a few RG steps all other operators 
are already in the FP, so the fiow line follows that single operator. According to 
Eq. [5] the scaling exponent of the gauge coupling for the Nf = 16 fiavor theory 
at the Banks-Zaks FP is small, 7^ = —0.01, it is almost a marginal operator. The 
scaling exponents of the other irrelevant operators are likely close to their engineering 
dimensions, starting at 7 ~ —2, so these operators will die out much faster than the 
gauge coupling. 
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The picture we expect in the 2-lattice matching is now clear. Since the gauge couphng 
is nearly marginal, the matching will follow its flow. Sf,(/5) is given by Eq. [7] and 
for a marginal operator Sb{(3) = near the FP. For an almost marginal operator 
A/5 = (3 — (3' can be either positive or negative, depending on whether the FP of the 
actual RG transformation is at smaller or larger gauge coupling. 

The evolution of the blocked operators can also signal the IRFP. As the RG flow 
approaches the IRFP all expectation values approach their IRFP value. On the other 
hand if the GFP controls the system the flow lines follow the RT and run into the 
trivial /5 = FP where all local expectation values vanish. This difference is one of 
the strongest signal for the existence of an IRFP in the MCRG method. 

• At flnite mass the RG transformation is dominated by the flow of the relevant mass 
operator. However the nearly marginal gauge coupling can still have a strong influence 
on the flow, the situation is more like matching 2 relevant operators than matching 
a single one. The easiest way to deal with this is to set the gauge coupling to its 
FP value (i.e. to the value that corresponds to the IRFP of the RG transformation 
used) and match in the mass only. This matching predicts the scaling dimension (or 
critical exponent) of the mass. While the IRFP of the RG transformation depends on 
the blocking parameter, the exponent itself is independent of both a and the gauge 
coupling (3. 

2. Numerical simulations and results 

I concentrate on the renormalization group properties of the model in the massless limit 
and show only preliminary results at flnite mass. I did 2-lattice matching on 16^ — > 8^ lattices 
using nHYP smeared staggered fermions [43], and as in the Nf = 4 case, the simulations 
were done with a small mass, m = 0.01 on the 16"^ and m = 0.02 on the 8^ lattices. I have 
collected 100-150 conflgurations on the larger volumes, ~ 200 on the smaller ones, separated 
by 10 molecular dynamics steps. 

On the 8^ lattices I covered the coupling range f3 G (2.4,8), trying to identify the bulk 
transition to the strong coupling phase. The data show no sign of a phase transition, 
though for P < 4.0 the condensate starts increasing slowly, suggesting the development of 



23 




0.0 



0.2 



0.4 



m 



FIG. 11: The dependence of the plaquette on the mass after 1 (red diamonds and sohd line) and 
2 (blue diamonds and solid line) blocking steps on 8^ volumes at (3 = 5.8. The bursts at m = 0.15 
show the plaquette after 2 and 3 blocking steps starting form 16^ volumes at the same P value. 
The dashed lines indicate matching in the mass. The blocking parameter is a = 0.75, close to the 
optimal value at (3 = 5.8. 

spontaneous chiral symmetry breaking. It is likely that the bulk transition exists only at 
very small, possibly vanishing, mass. This does not contradict the results of Ref. 471] where 
a strong first order bulk transition was observed with Nf = 16 fermions. As we learned 
from the finite temperature investigation with Nf = 4 flavors in Sect. IIVB 11 smearing the 
fermionic action can soften or wash away first order phase transitions, and that might be 
the situation here as well. 

Just like in the Nf = 4 case the mass dependence of the measured operators is weak for 
small m. Figure [11] shows the plaquette blocked 1 and 2 times on the 8^ configurations at 
/3 = 5.8 at different masses. The solid lines are simple spline extrapolations. Within errors 
there is no mass dependence up to about m < 0.05 for the plaquette. Other observables are 
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FIG. 12: The matching of the plaquette for Nf = 16 flavor. The individual data points are at 
(3 = 5.6 on 16^ lattices after 2 and 3 levels of blocking. They are compared to the 1 and 2 times 
blocked values as measured on the 8^ lattices (solid lines). The matching values A/? ~ indicate 
a nearly marginal flow. Note that the plaquette increases with rif, implying that the RG flow is to 
an IRFP and not the /3 = trivial FR 

similar, so for the MCRG matching the m = 0.01 — 0.02 data set can be considered critical 
in the mass. 

Figure [12] is the analogue of Figures [2] and [HI showing the matching of the plaquette. 
The data points are a.t (3 = 5.6 on 16^ lattices after 2 and 3 levels of blocking. They are 
compared to the 1 and 2 times blocked values as measured on the 8^ lattices. The block 
transformation is with parameter a = 0.78, optimal for (3 = 5.6 and matching is consistent 
for both blocking levels with Af3 = 0. 

From the matching value alone it is not possible to distinguish a marginally relevant 
flow (like at the GFP) from an almost marginal irrelevant flow (expected at the IRFP). 
Comparing Figure [12] to the Nj = 4 or pure gauge Figures [2] and [8] reveals an important 
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FIG. 13: Matching at jS = 5.4 (left panel) and f3 = 6.6 (right panel) from 16^ to 8^ lattices in the 
Nf = 16 flavor theory. 

difference. When the flow is governed by the GFP the flow lines follow the RT towards the 
trivial f3 = PP. In the — > oo limit all expectation values vanish. In Pigures [2] and M the 
plaquette indeed decreases with increasing blocking levels. On the other hand when the flow 
lines approach an IRPP all expectation values take the value at the PP. The plaquette in 
Pigure [12] increases with increasing blocking steps for f3 > 5.4, indicating that the flow lines 
are not running towards the (3 = PP. Eventually all expectation values should become 
independent of the blocking steps and (3. This trend is not obvious from the flgure for two 
main reasons. The flrst is that flnite volume effects are considerable on 2^ lattices, the second 
that with only 2 or 3 blocking steps one approaches the PP only if the RG transformation is 
optimized. Nevertheless the fact that with optimal blocking the blocked operators increase 
with rif, already signals that the flow runs toward an IRPP. 

Other operators show similar behavior at /5 > 5.6 but the results are quite different at 
stronger gauge couplings. Pigure [12] shows that the blocked plaquette values cross around 
f3 = 5.4 and they decreases with Ub for (3 < 5.4. The location of the crossing depends on 
the blocking parameter but for the Nf = 16 flavor model it never drops below (3 = 5.4. It 
appears that the flow is running towards the j3 = PP when (3 < 5A and to the IRPP when 
f3 > 5.4. 
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/3l6 


Q^opt 


Sb{s = 2) 


5.4 




none 


5.6 


0.78 


0.00(3) 


5.8 


0.78 


-0.12(4) 


6.2 


0.66 


0.09(5) 


6.6 


0.67 


-0.02(4) 



TABLE III: The parameters and MCRG results of the Nf = 16 flavor simulations. The second 
column lists the optimal blocking parameter a. 

0.5 I — I 1 1 1 1 1 1 1 1 1 — I 1 1 1 — 




0.0 0.1 0.2 0.3 

mi 

FIG. 14: Matched (mi, 7712) pairs in the A'^^ = 16 system. Red diamond: (3 = 5.8, blue bursts: 
P = 6.6. The Unear fit predicts ym = 1.02(7). 

The 2-lattice matching is also different below and above (3 = 5.4. Figure [I3] shows A/3 
as the function of the blocking parameter a for (3 = 5.4 and (3 = 6.6. This figure is the 
analogue of Figure HI As the left panel of Figure [13] shows the = 2(1) and 3(2) blocking 
levels get close but do not actually converge at /? = 5.4, there is no consistent matching. 
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(3 


"opt 


mi 


m2 


5.8 


0.75 


0.10 


0.17(3) 


5.8 


0.75 


0.15 


0.31(2) 


6.6 


0.68 


0.07 


0.16(3) 


6.6 


0.68 


0.10 


0.20(4) 


6.6 


0.68 


0.15 


0.33(4) 



TABLE IV: Matched mass pairs in the Nj = 16 system. 

The situation is similar, even more enhanced, at stronger couphngs. I have not observed 
this kind of behavior either with Nf = or 4, though I had simulations at even stronger 
couplings there. While it is possible that higher blocking levels would predict matching 
and the expectation values eventually approach their FP value, it is more likely that we 
see the effect of a nearby bulk transition and the strong coupling confining region beyond 
it. Apparently (3 < 5.4 is not governed by the IRFP. The situation is entirely different at 
f3 = 6.6 where predictions from the two different blocking levels converge around a = 0.675 
predicting Af3 = —0.022(44) (right panel of Figure [T3!) . Results are similar at other f3 > 5.4 
couplings, as summarized in Table IIIII The optimal block transformation predicts A/3 ^ 
for all coupling values, in agreement with the expectations, i.e. that the RG flows are 
governed by an almost marginal operator. Combining this with the observation that with 
optimal blocking parameter the expectation values of the blocked operators increase leads 
to the conclusion that the (3 G (5.6, 6.6) coupling range in the massless limit is governed by 
an IRFP. 

The mass is a relevant operator at this IRFP, therefore it should scale according to Eq. 
El under an s = 2 RG transformation. One can calculate the exponent ym by identifying 
matched [1711,7712) mass values. The gauge coupling is irrelevant, any (3 in the attractive 
basin of the IRFP could be chosen for this. Since the gauge coupling is nearly marginal it 
is best to use the same coupling at both mass values. For matching one can use the same 
operators as before or add others that are more sensitive to the ferminons. I believe that 
direct fermionic observables would allow more precise matching at smaller quark masses, 
but as Figure [TT] illustrates the gauge observables also work. In fact Figure fTTl already shows 
a matched mass pair. At f3 = 5.8 with a = 0.75 RG transformation the mi = 0.15 mass 
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on the 16^ configurations match m2 = 0.318(6) on the 8"^ configurations after 2(1) blocking 
steps, while the matching mass is 1712 = 0.297(10) after 3(2) blocking steps. These values 
are for the plaquette but other observables give similar values, predicting 1712 = 0.31(3) at 
the optimal a = 0.75 parameter. I have preliminary data at /5 = 5.8 and 6.6 at a couple of 
mass values. The matching (mi, 7712) pairs are listed in Table HVl The fit according to Eq. [S] 
predicts t/m = 1.02(7) as shown in Figure The scaling dimension of the mass is very close, 
within errors undistinguishable, form it engineering dimension. This is not unexpected as 
the IRFP is at weak coupling. In a recent publication = 1.5 was predicted for sextet 
fermions 491]. It would be interesting to study the sextet model, or the Nf = 12 model, 
where the t/m might be significantly different from 1. 



V. CONCLUSION 



Renormalization group methods have been designed to study the critical behavior of 
statistical systems. In this paper I have shown that they are equally suitable to study the 
renormalization group structure of quantum field theories. I have used a numerical Monte 
Carlo Renormalization Group method to calculate the step scaling function and critical 
exponent of SU(3) gauge theories with Nf = 0, 4 and 16 fiavors. I chose these fairly well 
understood systems as my goal was to test the method before using it in more relevant 
simulations. The paper is fairly pedagogical, explaining in detail the 2-lattice matching 
MCRG method. 

In the Nf = case I demonstrated that the bare step scaling function predicted by the 
2-lattice matching method is consistent with the more traditional Schrodinger functional 
results. It is also consistent with results obtained from the scaling of the tq parameter and 
the finite temperature phase transition even at strong gauge coupling, suggesting scaling, 
though not 2-loop perturbative scaling there. 

The Nf = 4 and 16 fiavor simulations were done using nHYP smeared staggered fermions. 
I chose nHYP smearing because its highly improved taste symmetry. In the Nf = 4 fiavor 
case I briefiy studied the finite temperature phase transition to develop a feel for the pa- 
rameters of the model. I calculated the step scaling function at vanishing quark mass and 
showed that the 2-lattice method works equally well in the fermionic system. 

The Nf = 16 fiavor model brings in new challenges as it is governed by an infrared fixed 
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point at finite gauge coupling. At this IRFP the scahng dimension of the gauge couphng is 
small, it is an almost marginal operator. Accordingly the gauge coupling runs very slowly. 
The measured step scaling function is consistent with an almost marginal operator, but on 
its own it cannot predict if the gauge coupling is relevant or irrelevant. I argued that the 
evolution of blocked operators signal if the RG flow is towards an IRFP or to the (3 = 
trivial FP. For the Nf = 16 model the flow clearly indicates an IRFP. 

Finally I presented preliminary measurements for the scaling dimension of the mass. I 
found Urn = 1.02(7), undistinguishable from the free field exponent. This is not surprising 
for the Nf = 16 theory. 
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